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We investigate the spin relaxation due to the random Rashba spin-orbit coupling in symmetric 
GaAs (110) quantum wells from the fully microscopic kinetic spin Bloch equation approach. All 
relevant scatterings, such as the electron-impurity, electron-longitudinal-optical-phonon, electron- 
acoustic-phonon, as well as electron-electron Coulomb scatterings are explicitly included. It is shown 
that our calculation reproduces the experimental data by Miiller et al. [Phys. Rev. Lett. 101, 
206601 (2008)] for a reasonable choice of parameter values. We also predict that the temperature 
dependence of spin relaxation time presents a peak in the case with low impurity density, which 
originates from the electron-electron Coulomb scattering. 

PACS numbers: 72.25.Rb, 71.10.-w, 71.70.Ej, 73.21. Fg 



Semiconductor spintronics has been an active field of 
research lately due to the potential application of spin- 
based devices^ - — Recent experiments show that the spin 
relaxation time (SRT) in (llO)-orientcd GaAs quantum 
wells (QWs) is extremely long, and thus the spin dy- 
namics in this system has attracted much attention both 
experimentally and theoretically.— ~— The physics under- 
lying this effect is the absence of the D'yakonov-Perel' 
(DP) mechanism,— which is the leading spin relaxation 
mechanism in n-type zinc-blende semiconductors. The 
DP mechanism is from the joint effects of the momen- 
tum scattering and the momentum-dependent effective 
magnetic field (inhomogenous broadening- ) induced by 
the Dresselhaus^I and the Rashba^ 2 - spin-orbit coupling 
(SOC). In symmetric GaAs (110) QWs with only the 
lowest subband occupied, the in-plane component of the 
spin-orbit field vanishes^ Therefore the DP mechanism 
cannot affect electrons with spin polarization along the 
growth direction and the SRT in this system is consid- 
erably larger than that in (100) QWs. In most of the 
previous works, the main reason limiting the SRT is at- 
tributed to the Bir-Aronov-Pikus mechanism,— which is 
from the exchange interaction between the electrons and 
the photo-generated holes. One of the exceptions is the 
spin noise spectroscopy measurement by Miiller et al.^ 
where the excitation of semiconductor is negligible and 
hence the Bir-Aronov-Pikus mechanism is avoided. They 
reported the longest SRT in this system which is about 
24 ns. Since the DP and Bir-Aronov-Pikus mechanisms 
are both absent, and the virtual intersubband spin-flip 
SOC induced spin relaxation is also ruled out due to the 
relatively high mobility of the samples, 16 the possible rea- 
son limiting the SRT is the DP mechanism due to the 
random Rashba SOC caused by the fluctuations of the 
donor density 25 > 26 

As shown by Sherman et aL ) 25 ' 26 even in symmetric 
QWs, the unavoidable fluctuations of the concentration 
of the dopant ions still lead to a random electric field 
along the growth direction, and hence a random Rashba 



SOC at each point of a QW. This random SOC provides 
an inhomogeneous broadening and induces the DP spin 
relaxation. The previous investigation s 25 ! 26 on the spin 
relaxation due to this mechanism are based on the single- 
particle theory, and thus the electron-electron Coulomb 
scattering, which has been shown to be very important 
for spin relaxation in two-dimensional system , 20 i 27 ~ — is 
missing. In this work, we apply the fully microscopic ki- 
netic spin Bloch equation (KSBE) approac h 20 ! 28 to inves- 
tigate the spin relaxation due to the random Rashba SOC 
in symmetric GaAs (110) QWs. Here all relevant scatter- 
ings, especially the electron-electron Coulomb scattering, 
are explicitly included. We will show that our calcula- 
tion is in good agreement with the experimental data by 
Miiller et al.J^- We also predict that the temperature de- 
pendence of the SRT presents a peak due to the electron- 
electron Coulomb scattering in the low-impurity-density 
case. 

We start our investigation from an n-type symmetri- 
cally modulation-doped GaAs (110) QW with the growth 
direction along the z axis. The in-plane coordinate axes 
are set as a; || [110] and y || [001] . A magnetic field B is ap- 
plied along the a;- axis. The well width is small enough 
so that only the lowest subband is occupied for the tem- 
perature and electron density we discuss. The envelope 
function of the relevant subband is calculated under the 
finite-well-depth assumption; 2 — The barrier layer is cho- 
sen to be Alo.39Gao.6iAs as the experiment where the 
barrier height is 319 meV^ By using the nonequilib- 
rium Green's function method^ 5 - the KSBEs can be con- 
structed as3 

dtPk = <9tPk|coh + <9 t /3 k | S cat, W 

where pk represents the electron single-particle density 
matrix, whose diagonal and off-diagonal elements de- 
scribe the electron distribution function /k<r and spin 
coherence pk respectively. The scattering term dtPk\coh 
consists of the electron-impurity, electron-longitudinal- 
optical (LO)-phonon, electron-acoustic (AC)-phonon 
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and electron-electron Coulomb scatterings, whose expres- 
sions are given in detail in Ref. 28, The coherent term 
can be written as (H = 1 throughout this paper) 



dtPk 



coh 



h(k)-(T + S H F(k), p k 



(2) 



in which [A, B] = AB — BA is the commutator. EnF(k) 
is the effective magnetic field from the Coulomb Hartrcc- 
Fock contribution. 20 h(k) represents the spin-orbit field 
composed of the DresselhausSi and Rashba— terms. The 
Dresselhaus spin-orbit field in (110) QWs reads 



h D (k)= 7D (0, 0, ^{kl-2k 2 y -{kl)) 



(3) 




Here 7d denotes the Dresselhaus SOC coefficient and (k 2 ) 
stands for the average of the operator — (d/dz) 2 over the 
electronic state of the lowest subband. The effective mag- 
netic field from the Rashba SOC can be written as 



FIG. 1: (Color online) SRTs from the experimental data in 
Ref. [l^ (•) and from the KSBEs (solid curve) vs. the applied 
magnetic field. Here T = 20 K, a = 16.8 nm, iV e = 1.8 x 
10 11 cm -2 and JV< = 0.01 N e . 



hR.(k) = a R (k y , -k x , 0) 



(4) 



with an representing the Rashba SOC coefficient. To in- 
corporate the effect of the random Rashba SOC, we as- 
sume that q;r, satisfies the Gaussian distribution function 
P(a-n) — ^ /2 1 A e~ a rt/ 2A " with Ar, being the standard 
deviation of the distribution function.— We divide the 
regime of an, from — a cu t to a cu t equidistantly by JVr 
odd nodes, where a cu t is the cutoff value of an regime. 
For each node of an, we obtain the temporal evolutions 
of the electron distribution function by numerically solv- 
ing the KBSEs. The SRT is obtained by the slope of the 
envelope of the coherently summed spin polarization 

S z = J2 A«rP(«r) ]T i[/kt(a R ) - /hK)] 

«r k 

with Aor = 2a cu t/(-/VR — 1). It is checked that a cu t = 
1.6 meV-A and Nn = 15 is sufficient for the convergence 
of our calculation. 

First, we compare the calculation via the KSBE ap- 
proach with the experimental data in Ref. 15, In Fig. [TJ 
we plot the SRT calculated from the KSBEs together 
with the experimental data as function of the applied 
magnetic field. Here temperature T = 20 K, well width 
a = 16.8 nm and electron density N e — 1.8 x 10 11 cm -2 
as indicated in the experiment^ 5 Impurity density Ni = 
0.01iV e is obtained by fitting the transport mobility given 
in Ref. [lit It is shown that our calculation agrees well 
with the experimental data for various magnetic field 
strength. The fitting gives Ar = 0.195 meV-A and 
7d = 5.0 eV-A 3 . As shown by Sherman^ 5 - the random- 
ness of the Rashba SOC is related to the effective distance 
between the QW conducting sheet and the dopant layer 
zq. The relationship can be written as 
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FIG. 2: (Color online) SRTs due to the DP mechanism in- 
duced by the random Rashba field vs. temperature T for 
B = (a) and 4 T (b) with impurity densities iV, = N e (•) 
and 0.01 N e (□). Here a = 16.8 nm and N e = 1.8xlO n cm" 2 . 
The corresponding Fermi temperature is 75 K. The red solid 
curves represent the results calculated with all relevant terms; 
the blue dashed curves are those without the electron-electron 
Coulomb scattering; the green dotted curves are those with- 
out the cubic Dresselhaus SOC. 
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where £ = 7 e-A 2 for GaAsJ^ From the fitting value of 
Ar, the corresponding effective distance z = 21 nm, 
which is consistent with the sample parameter given in 
Ref. [Tip! Our best fitted value of 7D = 5.0 eV-A 3 is 
fairly close to the previous fitting value in GaAs (100) 
QWs (5.7 eV-A 3 in Ref. and 8.6 eVA 3 in Ref. M) 
and the value reported in ah initio calculations (6.4, 
8.5 eV-A 3 )^ In fact, the measured values of 7d vary 
from 5.7 to 34.5 eV-A 3 ^^ This is because the Dressel- 
haus SOC coefficient is affected by the interface inversion 
asymmetry4i 

From Fig. [1] it is also seen that the SRT decreases with 
the magnetic field. The underlying physics is as follow- 
ing. As first shown by Wu and Kuwata-Gonokamii^ in 
the presence of an in-plane magnetic field, the in-plane 
and out-of-plane spins are mixed and the Dresselhaus 
spin-orbit field provides an inhomogeneous broadening 
which leads to the spin relaxation and dephasing. There- 
fore the SRT decreases with B. For the case with small 
magnetic field satisfying ujp <C (711 — 7 j_)/2, 42 the effec- 
tive spin relaxation rate can be written as^ 

7|| + 7J- 7J - 7± L 4faj| 

ic ^— — — v 1 -*^^ () 

where 7^ = T zi and 7|| = Tyy with z and y being the 
principle axes of the spin-relaxation-rate tensor T. This 
definition is due to the fact that the z-axis and y-axis 
are not the principle axes of T, when the Rashba and 
Dresselhaus SOCs are both present ^ 7 ' 18 It is noted that 
7_l differs from T zz which is the spin relaxation rate ne- 
glecting the Dresselhaus SOC. The physics of this effect 
is similar to that of the influence of the applied magnetic 
fields the Rashba SOC acts as an in-plane magnetic field, 
and mixes the in-plane and out-of-plane spins, thus the 
Dresselhaus SOC can affect the out-of-plane spin relax- 
ation. Our calculation gives the ratio 711 /"fj_ = 7.5, which 
is consistent with the huge spin dephasing anisotropy ob- 
served in experimentsi^&i^ 

In Fig. [3J we plot the SRT due to the DP mechanism 
induced by the random Rashba field as function of tem- 
perature for B = (a) and 4 T (b) with impurity den- 
sities Ni — N e and 0.01iV e . The other sample param- 
eters as well as Ar are chosen to be the same as the 
previous c&seM- We first concentrate on the case with- 
out magnetic field [Fig. Ufa)]. It is interesting to see 
that the SRT exhibits a peak around the Fermi tempera- 
ture Tp = Ep/kB in the low-impurity-density case. This 
peak originates from the electron-electron Coulomb scat- 
tering which dominates the momentum scattering. The 
Coulomb scattering increases with increasing T when 
T < Tp and decreases with T when T > T F ^2dl 
thus the peak appears around the Fermi temperature. 
It is also seen that the SRT increases with increasing 
T without the Coulomb scattering. This is because 
the electron-AC-phonon and electron-LO-phonon scat- 
terings both increase with temperature monotonically. 
From Fig. HJa), it is also found that the SRT decreases 



with increasing T monotonically in the high-impurity- 
density case. This is because the electron-impurity scat- 
tering dominates the momentum scattering in this case. 
Since the impurity scattering depends on temperature 
weakly, the temperature dependence of SRT is mainly 
determined by the inhomogeneous broadening from the 
SOC. With an increase of temperature, electrons are dis- 
tributed at higher momentum states. This leads to the 
increase of the inhomogeneous broadening and thus a de- 
crease of the SRT. 

Then we turn to the case with B = 4 T [Fig. gfb)]. 
For large magnetic field, the SRT is determined by 
7cff = 711 o 7J " ■ In consistence with the anisotropy fac- 
tor 7 ||/7_l = 7.5, it is seen that the SRT in this case is 
three to four times smaller than that in the case without 
magnetic field. The temperature dependence is similar 
to the case with B = 0: the SRT presents a peak around 
Tp in the low-impurity-density case, while decreases with 
increasing T monotonically in the high-impurity-density 
case. The only difference is that the SRT decreases faster 
at high temperature in the case with high magnetic field. 
This is because the contribution to the inhomogeneous 
broadening from the cubic Dresselhaus SOC becomes sig- 
nificant for large well width and high temperature.— By 
comparing the SRT without the cubic Dresselhaus SOC, 
it is found that the effect of the cubic term increases 
with T faster than the linear Dresselhaus term and the 
Rashba term (both are linear). This accelerates the de- 
crease of the SRT at high temperature. It is also seen 
that the SRT without the Coulomb scattering for Ni = 
first increases then decreases and again increases with 
increasing T. This effect is due to the competition of 
the increase in the inhomogeneous broadening and the 
increase in the scattering. When T < 50 K, the contri- 
bution from the linear SOC term is dominant and hence 
the increase of the inhomogeneous broadening with ris- 
ing temperature is slower than that of the electron-AC- 
phonon and electron-LO-phonon scatterings. This leads 
to the increase in the SRT. For temperature between 
50 K to 200 K, the contribution of the cubic SOC term 
is comparable to that from the linear ones. Thus the 
effect of the inhomogeneous broadening increases with 
increasing T faster than the electron-AC-phonon scat- 
tering but slower than the electron-LO-phonon scatter- 
ing. By further noticing that the electron-LO-phonon 
scattering surpasses the electron-AC-phonon scattering 
when T > 100 K [see the dashed curve with square in 
Fig. Ufa)], one can understand that the SRT decreases 
with increasing T when 50 K< T < 100 K but increases 
when T > 100 K. 

In conclusion, we have investigated the SRT due to 
the DP mechanism induced by the random Rashba field 
in symmetric GaAs (110) QWs via the fully microscopic 
KBSE approach, where all the relevant scatterings, espe- 
cially the electron-electron Coulomb scattering, are ex- 
plicitly included. We show that our calculation is in good 
agreement with the experimental data. We also find that 
the Coulomb scattering makes marked contribution to 
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the spin relaxation. It is predicted that the tempera- 
ture dependence of the SRT exhibits a peak in the low- 
impurity-density case, regardless of the applied magnetic 
field. This peak is from the nonmonotonic temperature 
dependence of the electron-electron Coulomb scattering. 
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